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Abstract. In this paper, we propose an algebraic approach to upgrade a projective 
reconstruction to a Euclidean one, and aim at computing the rectifying homog- 
^ ' raphy from a minimal number of 9 segments of known length. Constraints are 

derived from these segments which yield a set of polynomial equations that we 
solve by means of Grobner bases. We explain how a solver for such a system of 
ly-v ' equations can be constructed from simplified template data. Moreover, we present 

f^ , experiments that demonstrate that the given problem can be solved in this way. 



1 Introduction 

A projective reconstruction can be computed from image correspondences alone, with- 



> 

u 

c^ \ out any information about calibration or pose of the cameras . But additional knowledge 

about the scene is required to subsequently recover Euclidean structure |7 1. The prob- 
lem of computing such a homography that upgrades a projective reconstruction to a 
, metric one has been treated in many publications, e.g. in [[TSl lSll. 

JC I In our case, the additional knowledge about the scene comprises a set of segments 

Qs . with known lengths, i.e. pairs of points with known distances between them. This 

Q^ I problem has akeady been addressed in lITTI . where parallel projection was assumed 

^O ' and hence an affine reconstruction was supposed to be known. 

■T-H \ The projective case was considered in ifTSl . The authors introduced the concept 

f^ . of a Quadric of Segments (QoS), defined in a higher-dimensional space by the set of 

m ' segments of known length, that can be computed linearly. Euclidean structure can be 

recovered from the QoS in closed form, exploiting the fact that all spheres intersect the 
plane at infinity in the absolute conic Q- Although the linear computability is certainly 
advantageous, the high number of segments (54 in the 3D case) required to define the 
QoS constitutes a main drawback. Given that not many segments of known length are 
C^ . available in various conceivable settings, we therefore seek a solution that requires only 

as few segments as possible. 

The constraints on the upgrading homography derived from these segments consti- 
tute a system of non-linear algebraic equations. Polynomial systems occur in various 
computer vision problems and many of them have been solved by means of Grobner 
bases ||4l|5][r5|. But there is no easy, straightforward method to solve general polyno- 
mial systems efficiently and robustly. Instead, each particular problem usually requires 
the manual design of a suitable Grobner basis solver. However, an automatic generator 
of minimal problem solvers was presented in ifTOl . 

In this paper, we show that it is possible to compute the homography that upgrades 
a projective reconstruction to a Euclidean one algebraically from only 9 segments of 
known length. The following section states the problem in detail. Section [3] briefly 
introduces the notion of Grobner bases, on which our approach is based. We explain 
how a special solver for the problem at hand can be constructed and applied in sectionH) 
Finally, experimental results are presented in section |5] 
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2 Problem statement 

The image projection of a scene point Xi, represented by its homogeneous coordi- 
nates Q, is denoted as Xi oc PXi, i.e. an a^ G M exists, such that Xi = aiPXi. Let 
us assume that we have measured n image coordinates Xi qt Xi, i ^ N, and computed 
a projective reconstruction (P, {Xi}), where P 9(: P and Xi 9^ Xi, from those points, 
such that Xi oc PXi- Aiming at recovering Euchdean structure from this, we have 
to find a non-singular matrix H G M^^'' which upgrades Xi and P by a projective 
transformation, such that PXi oc PH^^HX^ and Xi ex HXi. 

In order to reduce the number of unknowns in H, we fix the reference frame by 
choosing three points Xi that determine the origin, the x-axis and the xy-plane. From 
now on, let us assume that all points Xi have already been mapped by a suitable 
similarity transform such that there exist three points (0, 0, 0, 1)^, {xi, 0, 0, 1)^ and 
{xj,yj, 0, 1)^, Xi,Xj, {jj eM.\ {0}, which determine the coordinate frame. 

Mapping the point of origin to itself, (0,0, 0,1)^ oc H(0,0,0,1)^, points on the 
X-axis to the x-axis, {xi, 0, 0, 1)^ oc H(a:i, 0, 0, 1)^, and points in the xy-plane to the 
xy-plane again, {xj , yj,0, 1)^ oc il{xj, ijj, 0, 1)^, the transformation matrix has to be 
of the form 



H 



/ hi 


h2 


h3 


\ 





h^ 


h5 











he 





\ hi-hg 


hj 


hs 


^9 / 



(1) 



Two constraints on H easily arise from that. First, the projection matrix has to be 
invertible and therefore H must fulfill 

j^ det(H) = hihihehg. (2) 

Secondly, we fix the scale of H by demanding that for one point Xi, I < i < n 

1 = H4Xi = X4i (3) 

with H4 denoting the 4-th row of H and X4i being the 4-th element of Xi. 

Furthermore, we assume that there are N pairs of points, {Xi, Yi) which represent 
segments of known lengths di, such that ||Xi — Y^ || = rfj for alii — 1 ... N. Replacing 
Xi by HXi, yields the following constraint on H: 

3 
fi{h) = = ^(HiX,H4f. - U^X^UiYif - (H4l,H4f,)'rf? (4) 

where H; denotes the ^-th row of H. 

Equation (JUl constitutes a homogeneous polynomial of degree 4 in 9 variables, with 
97 terms in the general case. Therefore, at least 9 segments {Xi, Yi) are required to 
obtain the 9 equations, which determine H. Introducing an additional variable hiQ, we 
can rewrite the inequality ^ as equality [2 1 

0=1 — hihihghghio, (5) 

as well as 

0=l-IliXi. (6) 

Now, the problem is to solve the system of 771 = iV + 2 algebraic equations, iV > 9, 
in 10 variables hi, ... , hio, 

= /i(h) = ... = /™(h). (7) 
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3 Grobner bases 

Systems of polynomial equations can be solved efficiently by means of Grobner bases ||6l . 
F denotes the set of m polynomials F = {/i (h) , . . . , f,n (h) | fi (h) G K[hi, . . . , /i„] } 
in n variables h — {hi, . . . , /i„) over a field K. The ideal / — (F) generated by F is 
the set of all polynomial linear combinations 

■Mh)q,{h)\q,{h) e K[hi, . . . ,K]\ . (8) 

A Grobner basis is a special set of generators with desirable algorithmic properties. 
In particular, a Grobner basis of an ideal / has the same set of solutions as /. But 
similar to a system of linear equations after Gaussian elimination, the solutions of 
I can be easily identified in the corresponding Grobner basis w.rt. a lexicographical 
monomial ordering |6|. 

Theoretically, the Grobner basis can be computed from any generating set of / 
by a method called Buchberger's algorithm [6|. The basic mechanism is to take each 
pair (/i(h), /j(h)) from F, /i(h) / /j(h), compute its 5-polynomial (see appendix), 
reduce it by F and add the remainder to F if it is not zero. This is done until the 
S'-polynomials of all pairs in F reduce to zero. 

This problem is known to be EXPSPACE-complete in general Q . Nevertheless, 
much better bounds can be found for many cases that actually occur in practice and 
several well-known methods exist to improve the basic algorithm. 

However, computing a Grobner basis straightaway from the set of equations in- 
troduced in section 12] with floating point arithmetic is not practicable for two reasons: 
One obstacle is that even with improved versions of Buchberger's algorithm many S'- 
polynomials are constructed in vain as they finally reduce to zero and hence do not con- 
tribute to the final basis, merely slowing the entire Grobner basis computation down. 
Another difficulty arises as a result of accumulating round- off errors in floating point 
arithmetic during repeated reductions of S'-polynomials. Except for extremely simple 
examples, these round-off errors make it impossible to decide whether a particular co- 
efficient very close to zero should be considered as zero causing the cancellation of the 
corresponding term or not. Thus a special solver that is adapted to the given problem 
has to be created. 



4 Solver design and application 

We build our solver by computing a Grobner basis of a template system of equations 
first. These polynomials are generated as explained in section |2] but originate from a 
simplified set of segments {Xi, Yi) with integer coordinates. Calculations are done in 
a finite field Zp. During the computation, we record which of the pairs of polynomials 
taken from the intermediate basis form S-polynomials that do not reduce to zero during 
the process. In so doing, we get a computation template that contains only those pairs 
of polynomials that actually contribute to the final Grobner basis. In principle, this step 
has to be performed only once. 

Afterwards, we can apply that solver to compute a Grobner basis from a general 
polynomial system in floating point arithmetic, taking only the previously recorded 
pairs of polynomials into account, which speeds up this procedure. To avoid the above 
mentioned problematic effects of round-off errors, we could have memorized when 
which integer coefficient during the template computation gets zero as proposed in 1 17] 
and proceed accordingly. But since storing and reading this information takes a con- 
siderable amount of memory and time, we favor the more practicable way of simply 
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processing the template data once again, simultaneously with the polynomial system 
that we want to solve. Whenever a coefficient of the template data (computed in Zp) 
becomes zero, the corresponding floating point coefficient is set to zero, too. 

We assume that the sequence of operations to construct the Grobner basis is basi- 
cally identical for different sets of polynomials, given that the equations in those sets 
contain the same monomials and differ only in their coefficients |16|. With this as- 
sumption, we rely on the fact that Buchberger's algorithm and its improved variants 
do not consider the values of non-zero coefficients for the choice of critical pairs, the 
detection of unnecessary pairs or the selection of reductors. 

The crucial point in this scheme is to find an appropriate template system that is 
simple enough to be feasible yet general enough to be used for the original problem. 
The difficulty is in the fact that having identical monomials in the template polynomials 
and in the original is required to achieve the same sequence of computation but does 
not necessarily lead to the desired result. 

To generate the template set, we simplified the general problem by using small 
integers coordinates for all segments \\Xi — Yi\\ = di, such that also di ^ Z\ {0}. 
That means each segment has to fulfill 

d^ = \\X,-Y,f=a^ + b^+c^, (9) 

where a,b,c,d E 'Z, d ^ 0, thus forming a Pythagorean quadruple if a, 6, c G Z \ {0} 
or a Pythagorean triple respectively if a = 0, 6 7^ and c ^ 0. The possibility to 
use such a special polynomial system as a template for the general problem is justified 
by the fact that every triple in M.^ has a sufficiently precise scaled representation as a 
Pythagorean triple in Z'^ [T4l. Hence, we are looking for a generic Pythagorean case 
which is feasible to compute and at the same time implementable for a wide range of 
practically occurring systems originating from real coefficients. 

5 Experiments 

5.1 Solver generation 

Although we aim to generate a solver to compute H from the minimal number of 
segments A^ = 9, we conducted experiments to assess the influence of the number of 
segments to the presented method. 

The template data set from which the solver is built consists of N pairs of points 
{Xi, Yi) that are generated from Pythagorean quadruples as outlined above, and a ho- 
mography H. More precisely, iV — 1 Pythagorean quadruples {ai,bi,Ci,di), where < 
ai,bi, Ci, di < 50, are selected randomly, such that (oi, bi,Ci,di) ^ {saj, sbj, scj, sdj) 
for any scaling factor s £ Z, for alH, j = 1 . . . A^ and i ^ j. Euclidean coordinates of 
Xi are chosen randomly as well, whereupon the points Yi are calculated, such that 

X,-Y,^ iXu,X2^,X3^,lV - {Yu , Y2^ , Y3, , 1^ = (a„6„c„0)^. (10) 

In general, Xi and Yi coordinates are non-zero integers, < \Xji\ , \Yji\ < 100, 
except for the first 2 segments 

Xi = (o,o,o,i)T, yi = (Yii, 0,0,1)^, 

X2 = (Xi2, X22, 0, 1)^, Y2 = (Y12, ^22, ^32, 1)^, 

which are required to determine the x-axis and the a;y-plane as explained in section |2] 
H is generated according to equation ([T), with random non-zero integer elements < 
hk < 20 for all A; = 1 ... 9. 
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The resulting A^ pairs of distorted points 

(H-iX„H-iy,) = (l„f,) (12) 

yield A^ + 2 polynomials according to equations (|4|i, (|5]l and ^. A Grobner basis is 
computed from this system of equations in a finite field Zp with p — 332251314113 as 
this proved to be a sufficiently large prime number in earlier experiments. In that way, 
solvers are built for 20 different template data sets for each of 10 distinct values of A^. 

As outlined in section [3] the intermediate basis F grows until the S'-polynomials 
of all pairs (/i (h) , fj (h) ) in F reduce to zero, usually producing a large final Grobner 
basis in this way. Therefore, its reduced Grobner basis [^ is computed afterwards, 
from which the solution can be easily obtained and which is unique for the considered 
ideal and a given monomial ordering. 

The reduced Grobner basis for the ideal generated by the considered polynomial 
system derived from N segments consists of the following 13 simple equations gi, that 
only differ in the coefficients Cy . 

51 = = Ci,i/li + Ci^2 92=0 = C2,l/l2 + C2,2 

33 = = C3_i/l3 + C3^2 34 = = C4_i/l7 + C4_2 

05 = = C5,i/l8 + C5^2 36 = = C6,i/l9 + C6,2 

97 = = C7,ih4 + cj^2h5 9s = = cs,ihl + cs,2 (13) 

gg = = C9,i/lg + Cg^2 5l0 = = Cio,ift.io + '^10,2 

511 = = cii^ih^he + 011,2/110 512 = = £124/15/110 + c^ahe 

513 = = Cisj/lg/llO + Ci3,2/l5 

These equations yield 4 solutions for H, varying in the signs of /i4, h^, and /ig, that 
maps all points Xi either to 



(14) 



and Yi analogously, depending on whether both /15 and /ig are positive, /15 < 
and /ig > 0, /15 > and /ig < or both are negative. Hence, a 180° rotation 
of all segments (Xl ,Y^ ) about the x-axis corresponds to the set of segments 
{(X^ , Kj )}, and mirroring these sets on the xz-plane yields the other two sets 
of segments {(X^+\ r/"+^)} and {{X^+~\y}+''>)}. 

In contrast to the explained basic Buchberger's algorithm, sets of pairs of polyno- 
mials {(/i(h), /j(h))} instead of single pairs are reduced simultaneously during so 
called multi -reduction steps in the modified version ||3] of the algorithm that we use in 
our experiments. 

The so computed reduced Grobner bases from various template data sets and dif- 
ferent values of N vary only in the coefficients c^ of equations ( fT3] l. Table [T] shows 
that the Grobner basis computation is faster for higher number of segments in terms 
of required multi -reduction steps and computing time. But more than 50 segments do 
not speed it up further Unsurprisingly, the final basis contains more equations if more 
multi-reduction steps were necessary to compute it. 

Next, the generated solvers are tested to find the minimally required precision for 
floating point arithmetic. More exactly, we apply each solver to its template data set, 
but compute its reduced Grobner basis in floating point arithmetic this time as explained 
in sectionlD If the absolute values of the difference between the so computed segment 



X^++^ 


— {Xii,X2i,X:ii,l) , X- 


= (Xu, 


-X2^,X3,,lV, 


x^+-^ 


^iXu,X2^,~X3„lV or X^—^ 


= {Xu, 


-X2^,~X3,,ly 
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N 


Basis 


Multi- 


Time 




size 


reductions 


in msec 


9 


500 


28 


5170 


10 


424 


24 


3828 


12 


364 


19 


759 


15 


206 


17 


140 


20 


105 


17 


11 


25 


95 


14 


10 


50 


67 


12 


5 


100 


67 


12 


7 


200 


67 


13 


20 


500 


67 


16 


57 



Table 1 : Solver generation: Size of the com- 
puted (not yet reduced) Grobner basis, re- 
quired multi-reduction steps and computation 
time for different numbers of segments A^. 



N 


Precision 


Time 




in bit 


(mimsec.msec) 


9 


1088 


2:49.23 


10 


512 


1:24.16 


12 


448 


0:12.50 


15 


384 


0:02.59 


20 


192 


0:00.11 


25 


256 


0:00.10 


50 


256 


0:00.05 


100 


256 


0:00.07 


200 


256 


0:00.20 


500 


256 


0:00.57 



Table 2: Solver testing: Required float- 
ing point precision and computation time 
to calculate the correct solution for differ- 
ent numbers of segments A^. 



length and the true segment length is smaller than 10^^, computed and true segment 
length are considered to be equal and the used precision is therefore sufficient. Our 
implementations make use of the GNU Multiple Precision Arithmetic Library f\} to 
handle high precision floating point arithmetic as well as large integers. 

The required precisions for different N are given in table |2] As expected, a higher 
precision is required for a lower number of segments, but using more than 50 segments 
does not decrease the necessary precision further. 



5.2 Application to exact data 

Now, we apply the solver to various floating point data sets. Each data set comprises A^ 
pairs {Xi, Yi), and a homography H. The Euclidean coordinates of the points Xi and 
Yi are in general randomly chosen within a cube of side lengths 1 0. All coordinates 
Xji,Yji G M are chosen such that lOX ji, lOYji G Z, i.e. only coordinates with at 
most one digit after the decimal point are considered. Furthermore, < Xji, Yji < 10 
for all j — 1 ... 3 and X4i = Yji = 1 for alH = 1 . . . A^. The true segment length 
di = \\X.i ~Yi\\ varies. 

H is chosen in a way that H^^ maps 3 vertices Vi of the cube to Vi+Ti oc H^^ V/, 
where T is a random vector of independent zero-mean Gaussian components with typ- 
ical deviation 1, and all hk G M. The system of polynomial equation is obtained from 
N pairs (Xi = li~'^Xi,Yi = H'^Yij according to equations Q, © and ^. For 
different values of N, each of the 20 solvers created in section ISTl is applied to 20 dif- 
ferent data sets, to compute the respective reduced Grobner bases, and then the sought 
homographies H'. Subsequently, the distorted points Xi and are YI are upgraded to 
X- = H'Xj and Y/ = H'l; and the segment lengths d'^ = ||X- - Y/\\ are computed. 

The error is calculated as the standard deviation of the differences between the 
computed and true segment lengths divided by the average segment length (T{di — 
d'^)/ ii{di). However, solvers may fail under certain conditions, such that no solution 
can be obtained. In table [3] the percentage of successful computed solutions and the 
mean errors of those solutions are summarized for different numbers of segments N . 

In our experiments, the solvers failed for instance on 2 data sets due to the chosen 
H^^ and a pair {Xi^Yi), that caused one or more coefficients in the resulting poly- 
nomial to be zero. Hence, a floating point exception occuiTed when division by zero 
was attempted during reduction by this polynomial. However, this problem occurred 
only in less than 0. 1 % of all polynomials that we generated from test data. In all cases. 
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N 


Success 


Mean 




rate 


error 


9 


100% 


7.5e-37 


10 


100% 


3.4e-45 


12 


100% 


4.1e-70 


15 


100% 


3.4e-82 


20 


100% 


2.3e-129 


25 


100% 


1.4e-137 


50 


95% 


2.3e-139 


100 


95% 


1.2e-141 


200 


85% 


4.5e-144 


500 


85% 


3.2e-149 



Table 3: Solver application to ex- 
act data: Success rate and mean 
errors for different numbers of 
segments TV. 



N 


Success 


Mean 




rate 


error 


25 


4% 


0.26 


50 


48% 


0.11 


100 


49% 


0.11 



Table 4: Solver application to 
noisy data, a — 0.001: Success 
rate and mean eiTors for different 
numbers of segments A*'. 





20 different data sets 


E 




x/xx/xx/x/,/x/x/x/x/x 


9 




X//XXXX////X/X/X/XXX 


9 




X/XX/XX////X/X/X/XXX 


9 




X//X/XX////X/X///X/X 


12 




X//XXXX////X/X/X/XXX 


9 




X//XXXX////X/X/X/XXX 


9 


e 


X//X/XX/X//XXX/X/XXX 


8 


> 

'o 

!Z1 


X//X/XX////X/X////XX 


12 


X/XXXXX////X/X/X/XX/ 


9 


0^ 


X/XX/XX////X/X/X/XXX 


9 


^ 


X/XXXXX////XXX///X/X 


9 


-a 


X//X/XX////XXX///XXX 


10 




■^•xxxxx////xx////xx/ 


11 




X/XXXXX////X/X/X/X/X 


9 




X/XX/XX////X/X/X/XXX 


9 




X//X/XX////XXX/X/X/X 


10 




X/XXXXX////X/X//XX/X 


9 




X//X/XX^//^X/X/XXX// 


11 




X//X/XX////X/X///XX/ 


12 



Table 5: Solver application to noisy data, A'^ = 50, 
a = 0.001: Success (/) and failures (X) for different ap- 
plying each of 19 different solvers to each of 20 distinct 
data sets. The last column displays the number of suc- 
cessful computations per solver. 



where the solver could be successfully applied to compute a solution, the upgraded 
segment lengths nearly equal the true length. 

5.3 Application to noisy data 

Finally, we investigated the effect of noise in the data. We used the same data sets 
as above, but considered only those solvers and datasets for which the computation 
did not fail in the previous experiment. A random vector of independent zero-mean 
Gaussian components was added to each Xi and Yi, and the polynomial system was 
derived from this noisy data. 

Our experiments revealed that even for a very small noise standard deviation, the 
computation completely failed in many cases. Here, failing means that though a re- 
duced Grobner basis was computed, no real valued solution could be obtained from 
that because the reduced basis contained equations like = ci h^ + C2 with positive 
coefficients ci and C2. Selected results are illustrated in tabled For larger numbers of 
segments the failure percentage as well as the mean error decreases, such that on some 
data sets none of the generated solver fails for TV > 50 as illustrated in table |5] 

6 Conclusion 



In this paper, we proposed an algebraic way to compute a homography to upgrade 
a preliminary projective reconstruction to an Euclidean one by means of constraints 
derived from a minimal number of segments with known lengths. We believe, this 
could be useful in environments were prior camera calibration is impracticable but a 
small number of distances between points is known. 
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We have shown that it is possible to solve this problem using only 9 segments and 
demonstrated how a corresponding solver can be constructed from simplified template 
data. However, our experiments revealed that the presented method is very likely to 
fail on noisy data. This has to be investigated more thoroughly in order to be able to 
apply our technique to real world data. 
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Appendix: Notation 

We use the notations term and monomial as they are explained in [6J, i.e. given a 
polynomial ring K[xi,X2, ■ ■ ■ ,a;„], a monomial is a product of the form 



with non-negative integer exponents ai, a2, . . . a„. A term then denotes the product 
ttax'^ of a monomial and a non-zero coefficient Oq G K- 

An S-polynomial of a pair of polynomials (/, g) is computed as 

5(/,5) = :r^(LT(/))-V-x^(LT(g)-ig 

where x'^ = LCM(LM(/), LM(5)) is the least common multiple of LM(/) and lM{g). 
LM(/) denotes the leading monomial and LT(/) the leading term of / w.rt. a mono- 
mial ordering. 

References 

[1] GMP: The GNU Multiple Precision Arithmetic Library. 

Ihttp: //gmplib. org 

[2] Thomas Becker and Volker Weispfenning. Grobner Bases: A Computational 
Approach to Commutative Algebra. Graduate Texts in Mathematics. Springer, 
New York, Heidelberg, 1993. 

[3] Michael Brickenstein. Slimgb: Grobner Bases with Slim Polynomials. Reports 
on Computer Algebra 35, Centre for Computer Algebra, University of Kaiser- 
slautem, 2005. 

[4] Martin Bujnak, Zuzana Kukelova, and Tomas Pajdla. 3D Reconstruction from 
Image Collections with a Single Known Focal Length. In Proceedings of the 
12th International Conference on Computer Vision, 2009. 

[5] Martin Byrod, Klas Josephson, and Kalle Astrom. Fast Optimal Three View 
Triangulation. In Proceedings of the 8th Asian Conference on Computer Vision, 
2007. 

[6] David A. Cox, John Little, and Donal O'Shea. Ideals, Varieties, and Algorithms: 
An Introduction to Computational Algebraic Geometry and Commutative Alge- 
bra. Undergraduate Texts in Mathematics. Springer, New York, Heidelberg, 
1992. 



OAGM/AAPR Workshop 2013 (aTXw:1304.1876D 



[7] Richard Hartley and Andrew Zisserman. Multiple View Geometry in Computer 
Vision. Cambridge University Press, Cambridge, 2006. 

[8] Anders Heyden and Kalle Astrom. Euclidean Reconstruction from Image Se- 
quences with Varying and Unknown Focal Length and Principal Point. In Pro- 
ceedings of the IEEE Conference on Computer Vision and Pattern Recognition, 
1997. 

[9] Klaus Kiihnle and Ernst W. Mayr Exponential Space Computation of Grobner 
Bases. In Proceedings of the 1996 International Symposium on Symbolic and 
Algebraic Computation, 1996. 

[10] Zuzana Kukelova, Martin Bujnak, and Tomas Pajdla. Automatic Generator of 
Minimal Problem Solvers. In Proceedings of the 10th European Conference on 
Computer Vision, 2008. 

[11] David Liebowitz and Stefan Carlsson. Uncalibrated Motion Capture Exploit- 
ing Articulated Structure Constraints. InternationalJournal of Computer Vision, 
51(3):171-187,2003. 

[12] Jean Ponce. On Computing Metric Upgrades of Projective Reconstructions under 
the Rectangular Pixel Assumption. In Marc Pollefeys, Luc Van Gool, Andrew 
Zisserman, and Andrew Fitzgibbon, editors. Second European Workshop on 3D 
Structure from Multiple Images of Large-Scale Environments, Lecture Notes in 
Computer Science, pages 52-61 . Springer Berlin Heidelberg, 2001. 

[13] Jose I. Ronda and Antonio Valdes. Euclidean Upgrading from Segment Lengths. 
International Journal of Computer Vision, 90(3);350-368, July 2010. 

[14] P. Shiu. The Shapes and Sizes of Pythagorean Triangles. The Mathematical 
Gazette, 67(439):33, March 1983. 

[15] H. Stewenius, D. Nister, F. Kahl, and F Schaffalitzky. A Minimal Solution for 
Relative Pose with Unknown Focal Length. In Proceedings of the IEEE Confer- 
ence on Computer Vision and Pattern Recognition, 2005. 

[16] Carlo Traverso. Grobner trace algorithms. In Proceedings of the 1988 Inter- 
national Symposium on Symbolic and Algebraic Computation, pages 125-138, 
1988. 

[17] Carlo Traverso and Alberto Zanoni. Numerical stability and stabilization of 
groebner basis computation. In Proceedings of the 2002 International Sympo- 
sium on Symbolic and Algebraic Computation, ISSAC '02, pages 262-269, New 
York, NY, USA, 2002. ACM. 



